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The paper discusses coherent control of charge and spin states of a biexciton system in a quantum 
dot via coherent two-photon transitions. Rabi oscillations between the ground state of a quantum 
dot and the biexciton state, as well as oscillations between the two single exciton states induced by 
laser pulses with different circular or linear polarizations are studied. The effect of phonon-induced 
decoherence on these processes is described. System properties and driving conditions that lead to 
optimal coherent control are identified. It is shown that proper optimization allows one to control 
the two-qubit biexciton system via two-photon transitions with a high fidelity. 
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I. INTRODUCTION 

Developing efficient methods for optical control of the 
quantum states of carriers confined in semiconductor 
nanostructures is not only important from a purely sci- 
entific point of view but also vital for emerging technolo- 
gies. Quantum dots (QDs) have been proposed as an 
efficient source of on-demand entangled photons 1 . This 
application was experimentally realized using optically 
driven QD systems 2 -. Demonstration of Rabi oscillations 
of the exciton occupation 3 and conditional control of the 
biexciton system (two confined excitons with opposite 
polarizations )^ demonstrate the feasibility of coherent 
optical control of charge (orbital) degrees of freedom. Al- 
though this suggests that implementing quantum bits on 
such charge states might be possible the relatively short 
lifetime of charge excitations considerably restricts the 
feasibility of this solution. Instead, spin degrees of free- 
dom, which are stable over much longer time scales, are 
believed to be much more promising . Here, again, opti- 
cal methods allow one to control the system on picosec- 
ond time scales, that is, much faster than by electron 
spin resonance or by electric gating using exchange inter- 
action. A whole range of theoretical proposals for such 
optical spin control schemes^^ was recently followed 
by an experimental demonstration 1 ^. It has also been 
shown that entanglement generating two-qubit gates can 
be performed optically on confined spin qubits via vari- 
ous all-optical scheme s 1 1 i 12 . 

While the atomic-like discrete properties of optical 
transitions in QDs allow one, in principle, to implement 
a wide variety of quantum optical schemes, the solid 
state nature of the system introduces various decoherence 
channels that cannot be ignored when designing optical 
control schemes. One of the most important sources of 
dephasing in coherently driven systems is the coupling 
to lattice degrees of freedom (phonons) 13 i 14 ' 15 which im- 
poses additional constraints on the implementation of op- 
tical control protocols 16 . Typically, due to the dynamical 
character of the lattice response and to the highly struc- 
tured nature of the lattice reservoir, the requirement of 
avoiding strong lattice response restricts the parameters 



of the control fields to a narrow range which can only be 
determined by careful modelling of the control procedure 
in the presence of decoherenc o 17 i 18 . 

In order to achieve fast control with high fidelity and 
to provide enough flexibility for optimizing against de- 
coherence processes, new control schemes are sought for. 
This paper is devoted to a specific class of such schemes: 
quantum control of a four-level biexciton system using 
coherent two-photon transitions. The direct motivation 
for the present study are recent experiments in which co- 
herent two-photon transitions between the ground state 
and the biexciton state have been demonstrated^ and 
full pulse-area-dependent two-photon Rabi oscillations 
between these two states were induced 2 ^. Here, a com- 
plete theory of coherent two-photon processes in a biex- 
citon system in a QD will be presented. Such processes 
provide an additional degree of control of the biexciton, 
which can be viewed as the simplest semiconductor two- 
qubit system. They can be used to create entanglement 
with a single laser pulse, which is a valuable alternative to 
the schemes based on sequential transitions or two-color 
control. From the spin-oriented point of view, some of 
the transitions to be described here consist in a simulta- 
neous flip of the electron and hole spin. Both the coher- 
ent evolution in an optically driven four-level system and 
the effect of carrier-phonon coupling will be discussed. It 
will be shown that QDs with various spectral properties 
(positive vs. negative biexciton shifts) and under various 
driving conditions (sign and value of the frequency de- 
tuning of the laser field) are optimal for different control 
schemes. 

The paper is organized as follows. Section |TT] intro- 
duces the model of the system under study. Next, in 
Sec. IIII1 the perfect system evolution (without dephas- 
ing) is discussed, first for linearly polarized laser pulses 
fSec. HII A"|) . then for circularly polarized ones (Sec. lIII 5|) . 
Section HVl outlines the numerical and analytical (pertur- 
bative) methods used for the study of the phonon impact 
on the system evolution. In Sec. [Vl the system evolution 
including phonon-induced dephasing is discussed for var- 
ious system configurations. Sec. I VII concludes the paper 
with a summary and final remarks. In addition, the influ- 
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FIG. 1: The energy level structure of a confined biexciton and 
the schematic presentation of the optical couplings between 
the states for the biexciton Rabi oscillations (a) and for the 
exciton spin flip (b). 

ence of the exchange splitting of the single exciton states 
is discussed in the Appendix. 

II. THE SYSTEM 

The system under consideration is composed of a single 
QD coupled to a laser beam which may have either linear 
or circular polarization. According to the selection rules, 
a circularly right-polarized beam can only induce tran- 
sitions from the ground state |g) to the right circularly 
polarized exciton state \<J+) and from the left circularly 
polarized exciton state |cr_) to the biexciton state |XX) 
(blue arrows labeled 'er + ' in Fig. [l}. The other two tran- 
sitions (red arrows labeled 'cr_' in Fig. [T]) are allowed 
for a left-polarized beam. A linearly polarized beam is 
a superposition of both circularly polarized components 
and couples both exciton states with fixed circular polar- 
izations (|o+), |<7_)) to the ground and biexciton states 
(|g), | XX)). In this way, the system is modeled by a four- 
level 'diamond' structure, which is a generalization of 'A 
systems' and 'V systems' studied in quantum optics. 

It is assumed that the frequency of the beam is de- 
tuned from all the single-photon transitions, as shown in 
Fig. [TJ For pulses of picosecond durations, the excited 
states of confined carriers are irrelevant and may be dis- 
regarded. In most of the discussion, we will also neglect 
the exchange interaction that couples the two circularly 
polarized exciton states and turns them into a weakly 
split linearly polarized double t 21 ' 22 ! 23 . The energy of this 
coupling is usually below 100 [ieV and can be reduced or 
even cancelled by applying external fields 24 - 25 or by spe- 
cial treatment (annealing) of the sample o 26 ' 27 . In most 
cases, the small fine structure splitting does not affect 
the system dynamics on the picosecond time scales rel- 
evant for the present discussion. However, it becomes 
important if the single exciton states are coupled by long 
(spectrally selective) pulses, so that the fine structure be- 
comes spectrally resolved. This effect is discussed in the 
Appendix. The effects of finite exciton lifetime (of the 
order of 1 ns) will also be disregarded. 



Thus, the system is described by the Hamiltonian 

H = Hxx + #las + H p h + -Hint. 

Here, the first term describes the four-level biexciton sys- 
tem, 

H xx = E{\a+)(a+\ + |a_X<r_|) + (2£ + £ B )|XX)(XX|, 

where E is the energy of the single exciton states and Eb 
is the biexciton shift (assumed negative in the binding 
case, as in Fig. [TJ. The second term accounts for the 
coupling between the carrier states and the laser beam. 
For a linear polarization it has the form 

Aas l} = /(*) cos(tti) [(|g) + |XX))(X| + H.C.] , 
while for a circular (say, <r+) polarization it reads 

^ = /(t)cos(nt) (\g)(a + \ + |0<XX| + H.c.) , 

where f(t) is the envelope of the laser pulse amplitude 
(which is assumed to be real), f2 is the laser frequency 
and we define 'linearly polarized' exciton states |X), |Y), 
related to the circularly polarized states (i.e., angular 
momentum eigenstates) by \a±) = (|X) ±i\Y))/y/2. The 
third term in Eq. ([TJ) is the free phonon Hamiltonian, 

where b\.,bk are creation and annihilation operators for 
a phonon with a wave vector fc and LOk is the correspond- 
ing frequency. The last term describes carrier-phonon 
interaction, 

-ffint = 

(|X)(X| + |Y)(Y| + 2|XX)(XX|) £ gl (b h + &L fc ) , 

h 

where the coupling constants gk — <?I fc account for the 
deformation potential coupling between the longitudinal 
acoustic phonons and a confined neutral exciton. This 
coupling dominates over the piezoelectric coupling to the 
acoustic phonon branches if the electron and hole wave 
functions overlap strongly 28 . On the other hand, cou- 
pling to optical phonons is not important in the present 
study since both the typical frequencies of the system 
evolution and the magnitudes of the detunings are much 
lower than the frequencies of the optical phonons. 

For the pulse envelope, a Gaussian shape will be as- 
sumed, 

/(t) = _^ e *UJ , 

V 27TT 

where tq is a parameter defining the pulse duration and 

»-f *^ 

J-oo n 
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is the pulse area. 

Assuming Gaussian wave functions and neglecting the 
small correction resulting from the different localization 
widths of electrons and holes, the coupling constants are 
given by 2 ^ 



» = ( ff --^)y^ e_, ' ( ^ )/ ^ /4 . 

where cr e ,h are the deformation potential constants for 
electrons and holes, p is the crystal density, v is the nor- 
malization volume for phonon modes, c is the speed of 
sound, and l,l z are the confinement lengths in the QD 
plane and along the growth direction. The values used 
in the calculations arc a c — = 9 eV, p = 5350 kg/m 3 , 
c = 5150 m/s, I — 4.5 nm, and l z = 2 nm. 

The Hamiltonian is transformed to the 'rotating frame' 
by the canonical transformation defined by the unitary 
operator 

jj _ e int(|<r + )(<7 + |+|o-_>(c7_|+2|XXXXX|)-i(A/ft)fl 

where A = hVl — E is the detuning between the laser 
frequency and the single exciton transition energy and I 
is the identity operator. The transformed Hamiltonian is 

H = UHU^ + ih^U 1 = H xx + H las + H ph + H int , 
where 

H xx = A|g)(g| + (£B-A)|XX)(XX|, (1) 
= /(t)cos(fii) 

x [(e- l0t |g) + e' at |XX)) (X| + H.c] , (2) 

<s C) = /(t)cos(fit) 

x[e~^(|g)(a + | + |a_>(XX|)+H.c.],(3) 

and the contributions Hph and Hi nt remain unchanged. 

Since the dynamics induced by the laser pulse is 
slow compared to the frequency of the optical field it 
is possible to treat the system in the rotating wave 
approximation^ -, as is commonly done in the description 
of picosecond dynamics induced by laser fields in the op- 
tical or near- infrared range (that is, with femtosecond 
oscillation periods). One neglects in Eqs. ^ and © 
the very quickly oscillating terms containing exp(±2ittt) 
which average to null over a typical time scale character- 
izing the evolution of the system state. This leads to the 
final formulas for the exciton-laser coupling Hamiltonians 

<r ) = ^/(t)[(|g) + |XX))(X|+H.c], (4) 

and 

<f C) = "/(*) (IS) Wl + IO<XX| + H.c.) . (5) 
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FIG. 2: The spectral branches representing the instantaneous 
eigenstates of the unperturbed Hamiltonian Hxx + H^ ln ^ as 
functions of the pulse amplitude / (the adiabatic parameter) 
for A = E B /2 (a) and for A = E B /3 (b). In both cases, 
E B < 0. 

III. UNPERTURBED EVOLUTION 

In this section, various schemes of two-photon coherent 
control of a biexciton system are discussed on a purely 
quantum-optical level, that is, without phonon-induccd 
dephasing. The optically induced transitions depend on 
the polarization of the laser pulse. The following two sub- 
sections are devoted to the evolution driven by linearly 
and circularly polarized laser pulses, respectively. 

A. Linear polarization 

Let us start with the unperturbed evolution of the sys- 
tem, generated by the Hamiltonian Hxx+H^' [Eqs. §Q 
and Consider the diagram of instantaneous eigen- 

states of the relevant three-level Hamiltonian (exclud- 
ing the decoupled state |Y)) as a function of the pulse 
amplitude /. First, let us focus on the special case of 
A = Eb/2 [Fig. Ha)]. Then 

H XX = |(|+)(+| + l-X-l), 

where |±) = (|g) ± |XX))/-y/2. The driving field [Eq. gj] 
does not couple the state |— ) to the other states, hence 
this state is invariant under H xx + . Another, non- 
trivial invariant subspace is spanned by the states |X) 
and |+). The instantaneous eigenvalues along the two 
spectral branches belonging to this subspace are 

A ± (t) = y{l±Vl + 2[/(*)/Ap}, (6) 

where A_ corresponds to the branch originating from the 
state |X) [upper branch, blue line in Fig.[2fa)]. It is clear 
that the two branches are separated from each other by 
at least |A|. 

If the amplitude / of the laser pulse changes slowly 
enough in time the evolution generated by H xx may be 
found using the adiabatic theore m 20 ' 31 . In such case, 
the states |X) and |+) undergo an adiabatic evolution, 
with f(t) playing the role of a slowly varying parameter. 
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At the time t\, after the pulse has been switched off, 
the initial state is restored with the additional dynamical 
phase 



a± 



dt\±{t), 



(7) 



where to is the initial time (before the pulse was switched 
on). 

Assume now that the system is initially prepared in 
the state 

v |X)+*|Y) 



V2 ■ 

The state |Y) is decoupled from the laser beam and does 
not evolve. As a result, the state \a+) undergoes the 
transformation 

e m -\X)+i\Y) 



V2 



e M ~ /2 icos- 



• • a - i 
I sm cr_ 

2 1 



The phase a_ may be arbitrarily large. Moreover, for 
a fixed pulse shape, it is a monotonous function of the 
pulse intensity. Thus, by varying the pulse amplitude, 
the exciton can be coherently rotated between the two 
polarization states |<r+) and |cr_). The resulting occupa- 
tion of the state |<t_), obtained from numerical solution 
of the quantum evolution equation, is shown in Fig. [21(a) 
(red solid line) as a function of the pulse area 9. 

In the same way, since the state |— ) evolves only triv- 
ially, if the system is initially in the state 



V2 

it will undergo the transformation 



<•'•" A,/| ;,) 1 fcos^|g) 



+ isin||XX) ) , 



where (3 = a+ + A(*i - t )/(2h). Note that A+ -> A 
when the laser pulse is switched off, so that (3 is in fact 
independent of the choice of the initial and final time. 
Thus, when the pulse amplitude is increased the system 
oscillates between the ground and biexciton states, much 
like in the usual pulse-area dependent Rabi oscillations 
between the ground and single exciton states, induced 
by a resonant circularly polarized bearrA The biexciton 
oscillations are plotted in Fig. [Stb) (red solid line) . Such 
oscillations were indeed observed in an experiment^. 

In spite of some qualitative similarity to the Rabi oscil- 
lations in a two-level system, described by the universal 
function sin 2 (9/2), one can clearly see essential differ- 
ences. The two-photon oscillations are not strictly pe- 
riodic, especially for weak pulses, when the transition 
probability develops very slowly. In fact, for weak pulses 
the occupation of the other state grows as 9 A ~ I 2 , where 
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FIG. 3: (a) Two-photon polarization flip: the occupation of 
the state |cr_) as a function of the pulse area for A = Eb/2 = 
—2 meV (red solid), A = Eb/3 (green dashed), and A = 
Eb/6 (blue dotted), (b) Two-photon Rabi oscillations^: the 
occupation of the biexciton state as a function of the pulse 
area for A = -Eb/2 = —2 meV (red solid), A = -1.9 meV 
(green dashed), and A = —1.7 meV (blue dotted). In both 
figures, Eb — —4 meV and to = 5 ps. 



/ is the pulse intensity, as expected for a two-photon pro- 
cess (see Ref . l20t) . Moreover, it is clear from Eq. ^ that 
the rotation angle is a nonlinear functional of the pulse 
envelope. Hence, contrary to the usual Rabi oscillations, 
no universal area theorem exists for the final occupations. 

In the general case, A ^ Eb/2, there are no in- 
variant states (apart from |Y)) and all three states 
|g),|X),|XX) give rise to three non-degenerate spectral 
branches [Fig. Hfb)]. Now, in the adiabatic limit, the 
state |X) evolves as previously, except for a different 
value of the corresponding adiabatic eigenvalue, which 
now cannot be given in a simple analytical form. Thus, 
the exciton spin flip effect will take place also in this 
case, as shown in Fig. [21a) (green dashed and blue dot- 
ted lines). 

On the other hand, since the states |g), |XX) now be- 
long to non-degenerate branches, an adiabatic evolution 
starting from any of these states will end up in the same 
state, up to an irrelevant global phase. Therefore, the 
two-photon Rabi oscillations are suppressed when the 
driving field is detuned from A = Eb/2, as shown in 
Fig- E^b) (green dashed and blue dotted lines). In fact, 
the symmetric detuning condition is equivalent to the 
two-photon resonance between the ground and biexciton 
states. If the two-photon resonance condition A = E/2\s 
not satisfied, the transition from the ground to biexciton 
state is forbidden and the corresponding Rabi oscillations 
are precluded. 

On the contrary, the polarization flip is always reso- 
nant and the detuning becomes a tunable parameter. As 
will be shown in the next section, this additional freedom 
of control may be useful for optimizing experimental pa- 
rameters against phonon-induced dephasing. 



B. Circular polarization 

In the case of a <r+ -polarized beam, the unperturbed 
system evolution is generated by the Hamiltonian Hxx + 



H 



(circ) 
las 



which has two invariant two-dimensional sub- 
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FIG. 4: (Color online) Two-photon rotation of the linear po- 
larization of an exciton state with a circularly polarized pulse: 
the occupation of the | Y) state as a function of the pulse area 
for Eb = — 4 meV and two values of A: = — 3 meV (a) and 
+2 meV (b). Red solid lines: to = 5 ps, green dashed lines: 
ro = 10 ps. The initial state is |X). 



spaces spanned, respectively, by the states |g), |<r+) and 
|XX),|cr_). In the adiabatic limit, each of these four 
states follows its own branch of instantaneous eigenstates 
and is restored after the laser is switched off, up to a 
phase factor. Like in the previously discussed cases, these 
phase factors can result in a nontrivial evolution if one 
starts with a system prepared in a superposition of the 
states. For dcfinitcncss, let us choose the initial state 



will be followed by the results of numerical simulations 
discussed in the next section. 

It is convenient to perform a unitary transformation to 
the phonon-dressed exciton states 3 ^, 



= E) 



(|X)(X| + \Y)(Y\)W + |XX)(XX|W 2 , 



where 



W = exp 



k HuJk 



(bk - bl k ) 



Note that this transformation conserves the occupations 
of the basis states. Working in the dressed basis physi- 
cally corresponds to the assumption that the initial state 
was prepared some time (a few ps) before the operation 
and it is therefore surrounded by the polaron-like coher- 
ent phonon field. Upon this transformation and expan- 
sion to the leading order in the coupling constants, the 
Hamiltonian may be written as H' = WifW + = H Q + V, 
where the the first component in the two cases of linear 
and circular beam polarization is, respectively, 



(lin) _ 

~~ 



A|g)(g| + H ph + {Eb - A)|XX)(XX| 
+%f(t) [(|g) + |XX))(X|+H.c] 



|X> 



\a+) + \a. 



V2 



The instantaneous eigenvalues corresponding to the two 
branches departing from |er + ) and |<r_) are 



A+ = A 



(8a) 



A = (Ey-A) 1 -^^*-^. (8b) 

The final state resulting from the adiabatic evolution is 
then 



cos ■ 



a.- , , olj 
|X) - sin — 



-|Y) 



where a± is defined by Eq. ([7]) but with X± defined in 
Eqs. (|5al b). Thus, the polarization of the exciton state 
can be rotated as long as a + — a_ ^ 0. This condi- 
tion is clearly not satisfied at the biexciton resonance, 
A = Eb/2, but, as long as the laser beam is detuned 
from this resonance, one expects polarization rotation to 
appear. This is confirmed by the numerical integration 
of the evolution equation, as shown in Fig. [5] Note that 
the final state is always linearly polarized, at an angle of 
(a + — a _) i '2 to the initial polarization. 



IV. PHONON EFFECTS: THE METHOD 

This section presents the general theoretical framework 
for calculating the effect of the phonon-induced decoher- 
ence on the two-photon transitions described above. This 



and 



H, 



(circ) 



- A|gXg|+if ph + (BB-A)|XX)(XX| 
+ f/W(|g}(^ + l + IO(XX|+H.c.) 



while the second term, describing the interaction with the 
phonon reservoir, can be written as V — ,$(W clrc 
where 



--^f(t)\-yx\ 



H.c. 



and 



s"" c) = -^/(*M»+I + k-XXX| + H.c, 



k 



Here. 



w = 1 



It. 



fjk 



(2n fe + 1) 



where nk is the phonon occupation number, accounts 
for the phonon-induced renormalization of the pulse am- 
plitude in the slow driving limi t 15 i 33 . We neglect the 
phonon-induced energy shifts which are very small for 
acoustic phonons. 

The evolution of the reduced density matrix of the 
biexciton subsystem is found by numerically solving the 
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time-convolutionless (TCL) evolution equation 3 ^ for the 
density matrix in the interaction picture with respect to 

p(t) = - / dr Tr ph [V(t), [V(r), p(t) ® Pph ]] , (9) 
Jo 

where V(t) is the interaction Hamiltonian V in the inter- 
action picture with respect to Hq, p p h is the phonon den- 
sity matrix at the thermal equilibrium, and Tr p h denotes 
the partial trace over the phonon degrees of freedom. 

The quantitative results presented in the following sec- 
tion will be obtained from a numerical solution to Eq. ((9]). 
However, much additional insight can be gained from a 
perturbative approximation to this equatio n 15 ' 35 ' and 
from the spectral interpretation it provide d 14 ' 37 . Let us 
note that, as long as the perturbation to the system 
evolution remains small, the reduced density matrix in 
the interaction picture differs little from its initial value. 
Hence, pit) on the right-hand side of Eq. ((9J may be re- 
placed by its initial value po = |V'o)(V'o|, where we assume 
that the initial system state is pure and represented by 
the state vector IV'o)- Upon integration, Eq. ([9]) then 
yields 

p{t) =po- f dr f dr' Tr ph [V(t), [F(t'), Po ® Pph]] • 
Jo Jo 

As a simple measure of the phonon-induced pertur- 
bation we will use the probability that the system in 
its actual state Uop(t)ul will be found in the desired, 
unperturbed state E/olV'o) (here, Uq is the unperturbed 
evolution generated by Ho). This is equal to T 2 = 
(ipo\p(t)\i>o), where T is the fidelity for the special case of 
a pure stated. It is convenient to write T 2 = 1 — 6, where 
5 is referred to as the error of a quantum control oper- 
ation. A physically transparent and meaningful formula 
for 5 is obtained by defining the frequency-dependent op- 
erator 

Y(u) = I [ drS{T)e luJT 
ft Jo 

and the phonon spectral density 

R(u) = f dt(R(t)R)e luJt , 

27T J-oo 

where S(t) and R(t) denote the operators in the interac- 
tion picture with respect to Hq. Then, one can write 

5= r (ii) 

where the spectral functions 

S l (u) = \& l \Y(cj)W )\ 2 (12) 

can be identified with different decoherence channels and 
the summation over / runs through all states orthog- 
onal to |-0o) (see Refs. mH for details). 




FIG. 5: (Color online) The two-photon Rabi oscillations in 
the presence of phonon-induced perturbation, (a) The de- 
pendence of the damping on temperature for tq — 5 ps: 
T = (red solid), T = 10 K (green dashed), and T = 40 K 
(blue dotted). The grey dash-dotted line shows the joint oc- 
cupation of the single exciton states due to real transitions 
at T = 40 K. (b) The dependence on the pulse duration at 
T = 4 K: to = 5 ps (red solid), to = 20 ps (green dashed) and 
to = 40 ps (blue dotted). 

V. PHONON EFFECTS: RESULTS 

This section is devoted to the system evolution for dif- 
ferent structures of the biexciton spectrum and under 
various driving conditions. In the first two subsections 
we will study the phonon impact on the biexciton Rabi 
oscillations and exciton spin flip induced by a linearly 
polarized laser field in a system with a negative biexci- 
ton shift (bound biexciton). The third subsection con- 
tains a discussion of the same processes in a system with 
antibound biexcitons (a positive biexciton shift) and the 
last one deals with phonon-induced dephasing in the two- 
photon rotation of the linear exciton polarization induced 
by a circularly polarized laser light. 

A. Biexciton Rabi oscillations for bound biexcitons 

Let us start the discussion with the case of biexciton 
oscillations in a system with a negative biexciton shift, 
Eb < (bound biexciton). In this case, the laser beam is 
detuned down from the transition to single exciton states, 
as shown in Fig.QJa). Therefore, one can expect that real 
phonon-assisted transitions to these states will be sup- 
pressed at low enough temperatures and the decoherence 
will be dominated by pure dephasing, like in the two-level 
cas o 14 ' 15 . Indeed, the oscillations presented in Fig. [5ja) 
show symmetric damping up to temperatures of several 
Kelvins. Only when A ~ fcrjT, the phonon-assisted tran- 
sition up to the single exciton state |X) is possible. As a 
result, the damping becomes much stronger and biased 
towards lower biexciton occupations, which is accompa- 
nied by a growing single exciton component, as shown by 
the grey dash-dotted line in Fig. 03a), corresponding to 
T = 40 K. 

The dynamical pure dephasing effect is related to the 
lattice relaxation after a non-adiabatic (with respect to 
the lattice response times) change of the charge stated 
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(a) ha>[meW] (b) hco[meV] 



FIG. 6: (Color online) The spectral density of the phonon 
reservoir at T = 4 K (thick grey lines on both plots, arbitrarily 
scaled) and the spectral functions Si (a) and S2 (b). Red 
solid lines: Eb = —4 meV, to = 5 ps; green dashed lines: 
Eb ~ —4 meV, To = 10 ps; blue dotted lines: Eb = —6 meV, 
to = 5 ps. Here always A = Eb/2. 

which correlates the exciton system with the phonon 
reservoir, generating a kind of "which path" information 
in the macroscopic environment^. Therefore, this kind 
of decoherence decreases when the optically driven sys- 
tem evolution becomes slower>2§,. This behavior is clearly 
visible in Fig. dfb). 

The properties of the two contributions to dephas- 
ing (pure dephasing and real transitions) may be conve- 
niently studied using the approximate perturbative for- 
mula (jTUjl . Fig. [5] shows the phonon spectral density 
R(oj)/ui 2 at T — 4 K and the two spectral functions 
Si(w), with |V>i) = I XX) and \ip 2 ) = |X), for two dif- 
ferent pulse durations. In all these plots, the nominal 
pulse area was tuned to achieve a complete Rabi flop [the 
first maximum in Fig. [3fb)]. One can see that the first 
function is always positioned around the zero frequency 
and its area (hence its overlap with R(lo)/uj 2 ) decreases 
for slower driving. Therefore, this spectral function may 
be identified with the pure dephasing effect which de- 
creases for slower driving. The shape of this function does 
not change considerably when the value of the biexciton 
shift Eb is altered. On the other hand, the second spec- 
tral function is located in the negative frequency area, 
around hu> ~ A. This allows us to relate this function 
to the real phonon-assisted transitions to the single ex- 
citon state |X). Since the phonon spectral density at 
large negative frequencies is small at low temperatures 
this decoherence channel is of minor importance in a sys- 
tem with a large enough biexciton shift and at sufficiently 
low temperatures. 

One notes in Fig. [6fb) that the shape of this spec- 
tral feature depends both on A and on tq. In order to 
study this real transition process more quantitatively let 
us note that for a localized spectral function S 2 (uj), the 
corresponding contribution to the error may be approxi- 
mated by 8 2 ~ hA 2 R(A)/A 2 , where 

/oo 
duiS 2 (u). 
-00 

The values of A 2 as a function of the pulse duration to 
are plotted in Fig.[7Ja). Again, for each t , the pulse area 
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FIG. 7: (Color online) (a) The value of A 2 /\A\ for E B = 
—2 meV (red solid), Eb = —6 meV (green dashed), and Eb = 
-8 meV (blue dotted). Here always A = Eb/2. (b) The 
first maximum of the biexciton occupation for A = Eb/2 = 
—2 meV, T — 4 K, and pulse durations to as indicated. The 
grey line is added to guide the eye through the maxima of the 
plotted curves. 



is adjusted in order to achieve a complete Rabi flop [the 
first maximum of the biexciton occupation in Fig. [3jb)] , 
so that all the points correspond to the same transfor- 
mation of the system state. The overall probability of 
a transition to the single exciton state depends not only 
on the duration of the process but also on the occupa- 
tion of the single exciton state during the evolution. The 
latter decreases for longer and weaker pulses. It turns 
out that the two effects compensate each other almost 
exactly, leading to a nearly constant value of A 2 as to 
is varied. This means that the error resulting from the 
real transition process is constant and, consequently, the 
biexciton occupation should be achieved with a constant 
(independent of t ) accuracy as soon as the evolution is 
slow enough for the pure dephasing to be negligible. This 
is confirmed by the numerical solution of the TCL equa- 
tion ([5]), presented in Fig. [7£b). An additional property, 
clearly seen in Fig. E[a), is that the value of A 2 is nearly 
exactly proportional to the detuning A, with just a small 
variation (a few per cent) for very short pulses. 



B. Exciton spin flip 

The spectral relations are different in the case of a two- 
photon spin flip. Now, for a negative biexciton energy, 
the laser frequency is always detuned up from at least one 
of the transitions [Fig. H(b)]. For the detuning as in the 
previous case, both transitions involve phonon emission 
and are therefore expected to cause much decoherence 
even at low temperatures. Indeed, the oscillations of the 
exciton spin orientation shown in Fig.[8ja) are damped in 
an asymmetric way (biased towards lower occupations), 
which is a signature of occupation leakage out of the two- 
dimensional single exciton subspace. Correspondingly, 
single exciton occupation becomes non-zero even at zero 
temperature. Moreover, the degree of decoherence does 
not decrease with a growing pulse duration [Fig. (SJb)] , 
again showing that the Markovian real transition pro- 
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FIG. 8: (Color online) The exciton spin oscillations in the 
presence of phonon-induced perturbation for Eb = 2A = 
—4 meV. (a) The dependence of the damping on tempera- 
ture for ro = 5 ps: T — (red solid), T — 10 K (green 
dashed), and T — 40 K (blue dotted). The grey dash-dotted 
line shows the joint occupation of the single exciton states due 
to real transitions at T = K. (b) The dependence on the 
pulse duration at T = 4 K: to = 5 ps (red solid), to = 20 ps 
(green dashed) and ro = 40 ps (blue dotted) . 




FIG. 9: (Color online) The spectral density of the phonon 
reservoir at T = 4 K (thick grey lines on both plots, arbitrarily 
scaled) and the spectral functions Si (a) and S2,S3 (b), for 
A = —2 meV (solid red), A = —3 meV (green dashed), and 
A = +2 meV (blue dotted). In all the cases, Eb = —4 meV 
and ro = 5 ps. 



cesses are a non-negligible factor. 

Interestingly, comparing Fig. [5] with Fig. [5] one can see 
that the overall dephasing of the two-photon spin flip is 
weaker than that of the biexciton oscillations, except for 
the first maximum at low temperatures. The reason for 
this is that the carrier-phonon interaction is insensitive to 
the spin orientation, hence there is no phonon-response 
to the transition between the two single exciton states. 
There could be some perturbation resulting from the oc- 
cupation of the other two states during the evolution but 
in the case of symmetric detuning, A = -Eb/2, both these 
occupations are equal and the average charge confined in 
the dot remains constant, which precludes any phonon re- 
sponse that might lead to pure dephasing. The spectral 
function related to the pure dephasing process is again 
Si(u), corresponding to |i/h) = |o+) in Eq. (fT2|). In the 
present case, this function remains null. This changes 
slightly if the detunings from the ground and biexciton 
states are not equal and, consequently, there is a charge 
variation during the evolution. However, as shown in 
Fig. [5J even in this case the spectral function Si(u>) is 



FIG. 10: (Color online) (a) The diagram of energy levels and 
laser frequencies in the case of A > 0, Eb < 0. (b) The two- 
photon exciton spin oscillations in the presence of phonon- 
induced perturbation for A > 0, Eb < for ro = 5 ps: T = 
(red solid), T = 10 K (green dashed), and T = 40 K (blue 
dotted). 

many orders of magnitude smaller than in the case of 
biexciton oscillations [compare Fig. [9ja) with Fig. [6ja)] . 

Phonon-assisted real transitions to the ground and 
biexciton states are now described by two spectral func- 
tions £2(0;) and £3(0;), corresponding to \1jj2) = |g) and 
|^3) = I XX) in Eq. (fl"2j) . respectively. As can be seen in 
Fig. EHb), these spectral features shift according to the 
detuning between the laser frequency and the relevant 
states. 

The fact that the spin flip does not require the two- 
photon resonance allows one to use the detuning A as a 
free parameter to optimize the control conditions against 
the phonon-induced real transitions. The most favorable 
choice is to set A > or A < Eb (still assuming Eb < 0). 
Then, one of the transitions is detuned down from res- 
onance and involves phonon absorption which makes it 
ineffective at low enough temperatures [see Fig. HOT a)]. 
The other transition is detuned up but the detuning can 
be chosen large enough to be brought beyond the cut- 
off of the phonon spectral density, as shown in Fig. [Hb) 
(dotted blue lines). The simulations based on the TCL 
equation ([9|), presented in Fig. fTOT b). show that spin flip- 
ping with fidelities of the order of 10 -3 is possible under 
such conditions. Note that this is achieved with a rela- 
tively short pulse so that the effect of the radiative decay 
during the operation is also small (of the order of 1%). 

C. Antibound biexciton 

Although the unperturbed evolution for Eb > (an- 
tibound biexciton) does not differ from that discussed 
above, this is no more true for the actual system kinet- 
ics in the presence of phonons. The two-photon reso- 
nance condition A = Eb/2 that must be satisfied for the 
two-photon Rabi oscillations to occur now means that 
the laser beam is detuned up from the single exciton 
transition. Unless this detuning is large enough (which 
would require large positive biexciton shifts) phonon- 
assisted transitions to the single exciton states may now 
take place, involving an emission of a phonon and there- 
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FIG. 11: (Color online) Dephasing of two-photon processes 
in a system with Eb > for to — 5 ps: (a) two-photon Rabi 
oscillations, starting from the ground state; (b) two-photon 
spin flip, starting from the \a+) state. Red solid: T = 4 K, 
green dashed: T = 10 K, blue dotted: T = 40 K. 



fore contributing to dephasing even at low temperatures. 
This can be seen in Fig. HiT a). where the results of nu- 
merical simulations are plotted at a few values of tem- 
perature. 

On the contrary, the positive biexciton shift is favor- 
able for the spin flip process since now both the tran- 
sitions to the ground and biexciton state are detuned 
down from the resonance. Hence, phonon-assisted tran- 
sitions may be expected to contribute considerably only 
at high enough temperatures. This is indeed confirmed 
by the numerical solution to the TCL equation shown in 
Fig. [TTT b). where it can be clearly seen that considerable 
dephasing occurs only when the temperature increases to 
a few tens of Kelvins. 

One should note that the strong phonon-induced de- 
phasing that affects the exciton spin flip process in the 
Eb < case could be avoided by detuning the laser beam 
strongly above the resonance, as discussed in Sec. IV Bl 
However, strong dephasing of the biexciton Rabi oscilla- 
tions in the present case cannot be dealt with in the same 
way because of the resonance condition required in this 
process. 



D. Exciton polarization rotation 

Let us now discuss the possibility of rotating the linear 
polarization of an exciton state via a two-photon process 
induced by a circularly polarized light. The energy di- 
agram for this process is the same as that relevant to 
the spin flip and is shown in Fig. [IJb). As discussed 
in Sec. IIII Bl performing this two-photon transfer re- 
quires detuning off the two-photon biexciton resonance, 
i.e., A 7^ .Eb/2. For a system with a bound biexciton 
(Eb < 0), this means that the detuning between the laser 
frequency and the transition to either the ground state 
or the biexciton state decreases and enters the range of 
high spectral density of the phonon reservoir. This re- 
sults in an increased contribution from phonon-assisted 
transitions and, therefore, in considerable damping of 
the polarization oscillations even at low temperatures, 



FIG. 12: (Color online) Dephasing of two-photon polarization 
rotation for to = 5 ps. (a) Eb ~ —4 meV, red solid line: 
A = -3 meV, T = 4 K; green dashed line: A = 2 meV, 
T = 4 K; blue dotted line: A = 2 meV, T = 40 K. (b) 
Eb = 4 meV, A = 3 meV, red solid line: T — A K; green 
dashed line: T = 10 K; blue dotted line: T = 40 K. 

as shown in Fig. [T2Ta). Efficient polarization rotation 
is only possible for large positive detunings [see the dia- 
gram in Fig. 110( a)], for the same reasons as discussed in 
Sec. IV Bl Again, the quality of oscillations is reduced as 
soon as the temperature becomes comparable with the 
detuning from the ground state. 

In a system with a positive biexciton shift, Eb > 0, 
phonon-induced transitions would require a phonon ab- 
sorption and are therefore suppressed at low tempera- 
tures, as shown in Fig.[T2Tb) (red solid and green dashed 
lines). It should be noted, however, that the system 
evolution is relatively slow in this regime of operations 
since the detunings are close to the two-photon reso- 
nance, when the effect of the laser field completely van- 
ishes (see Sec. IIII Bp . Therefore, large pulse intensities 
are needed to achieve the desired transformation of the 
system state. At higher temperatures [blue dotted line 
in Fig. 112( b)] the phonon absorption processes contribute 
considerably to the dephasing of the polarization rotation 
process. 

VI. CONCLUSION 

Two-photon transitions driven by picosecond laser 
pulses detuned from both exciton and biexciton reso- 
nances in a quantum dot open the way to a rich variety 
of control schemes in a biexciton system. Using a linearly 
polarized laser field one can induce Rabi oscillations di- 
rectly between the ground and biexciton state and co- 
herently flip the circular polarization of an exciton state. 
The former involves a two-photon absorption or emis- 
sion and requires exact two-photon resonance. The lat- 
ter can be interpreted as 'photon exchange' (absorption- 
emission) and can be induced in a broad range of detun- 
ings. With a circularly polarized field, one can induce a 
rotation of the linear polarization of an exciton. This pro- 
cess takes place only off the two-photon resonance. The 
physics beyond these two-photon processes is highlighted 
by invoking the quantum-mechanical adiabatic theorem. 
If the individual excitons are treated as separate quantum 
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subsystems (or qubits) then these two-photon control op- 
erations make it possible to prepare entangled states of 
these subsystems with a single laser pulse. 

Already these results, restricted to the ideal evolution, 
confirm once more that the biexciton system is an inter- 
esting object of investigation from the point of view of 
quantum optics. They show that couplings between all 
four states forming the 'diamond' structure of the biex- 
citon levels is essential for the system dynamics even if 
large detuning precludes real transitions between these 
levels. This suggests that models used for three-level 
atomic systems ('V-systems— ) may be insufficient for a 
description of QD systems. 

In a real semiconductor system, the evolution is 
strongly affected by phonon-induced dephasing. In ad- 
dition to the pure dephasing process appearing in an 
optically driven two-level exciton system, here we deal 
also with real phonon-assisted transitions to the states 
which, in the ideal case, should remain unoccupied. The 
effect of all the dephasing channels depends on the con- 
trol scheme under consideration and on the energy level 
structure of the biexciton system. For negative biexciton 
shifts (bound biexcitons), the phonon-induced decoher- 
ence is relatively weak for two-photon Rabi oscillations 
and becomes stronger for the polarization flipping pro- 
cess. However, the additional freedom of detuning allows 
one to optimize the driving conditions against these de- 
coherence effects also in the latter case so that spin flip 
fidelities as high as 10 -3 are achievable. QDs with a pos- 
itive biexciton shift are less suitable for these optical con- 
trol schemes. In conrast, they are favorable for the two- 
photon rotation of the linear polarization: in this case, 
optimization of the control conditions against phonon- 
induced dephasing can be much more efficient in a QD 
with an antibound biexciton than in a QD with a bound 
biexciton. 

The study of phonon-induced dephasing shows, in ad- 
dition, that a two-photon transition between the two- 
single exciton states can be performed via states with 
constant average occupation of the QD. This eliminates 
the phonon response to the charge evolution and allows 
one to avoid the resulting dephasing. With this respect, 
the single-pulse two-photon control has a clear advantage 
over sequential transitions. 

Out of the various processes discussed here, one 
(the biexciton oscillations) has already been observed 
experimentally—. The good agreement between the de- 
scription based on the adiabatic theorem and the experi- 
mental results shows that the theory captures the essen- 
tials of the quantum evolution under actual laboratory 
conditions. The other processes take place under sim- 
ilar experimental conditions and therefore also seem to 
be experimentally feasible. 




FIG. 13: (Color online) Rotation of the linear polarization of 
an exciton state with a circularly polarized pulse with the fine 
structure splitting of S — 30 /ieV (a) and S — 100 fieV (b): the 
occupation of the |Y) state as a function of the pulse area for 
Eb — —4 meV and A = — 3 meV; pulse durations 10 ps (red 
solid lines), 20 ps (green dashed lines) and 40 ps (blue dotted 
lines). 
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APPENDIX A: FINE STRUCTURE SPLITTING 

In almost all QDs, electron-hole exchange interaction 
splits the single exciton states into a linearly polarized 
doublet of states |X) and |Y). The magnitude 5{ s of 
this fine structure splitting can vary from several /xeV 
to about 100 /ieV. This splitting has no essential effect 
for the biexciton Rabi oscillations, while speaking of the 
exciton polarization flip makes sense only when the cir- 
cularly polarized states can be considered well-defined, 
that is, when the total duration of the experiment does 
not exceed H/6f s , which is of the order of 10 to 100 ps. 

In this Appendix the effect of the fine structure split- 
ting on the linear polarization rotation is discussed. The 
two linearly polarized states are well defined in the pres- 
ence of this splitting. Due to the requirement for energy 
conservation, the fidelity of this type of coherent control 
depends on the relation between the pulse length and the 
magnitude of the splitting. 

The additional contribution to the Hamiltonian, de- 
scribing the fine structure splitting, is 

if fB = ^(|X)(X|-|Y)(Y|). 

The results of simulations of the system dynamics (with- 
out dephasing) are presented in Fig. [13] for two values 
of the fine structure splitting 5f s . If the pulse is long 
and, therefore, spectrally narrow, the transition between 
the two non-degenerate states becomes forbidden by en- 
ergy conservation. As can be seen, for a typical value 
of 8f s — 30 fieV, the effect of lifting the degeneracy be- 
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comes important only for pulse durations of a few tens of picoseconds. 
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